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Abstract 

The problem of wave packet tunneling in a potential V{x) = (mu;^/2) (x^ — 5x'^) with 
I' > 2 is considered in the case when the barrier height is much greater than huj and the 
difference between the average energy of the packet and the oscillator ground state energy 
huj/2 is sufficiently small. The universal Poisson distribution of the partial tunneling rates 
from the oscillator energy levels is discovered. The explicit expressions for the tunneling 
rates of different types of packets (coherent, squeezed, even/odd, thermal, etc.) are given 
in terms of the exponential and modified Bessel functions. The tunneling rates turn out 
very sensitive to the energy distributions in the packets, and they may exceed significantly 
the tunneling rate from the energy state with the same average number of quanta. 
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1 Introduction 



Usually, the problem of quantum tunneling through potential barriers was considered under the 
assumption that the initial state possessed a well defined energy, i.e., for the quasistationary 
states. The propagation of the Gaussian wave packets through a rectangular barrier was studied, 
e.g., in [0. Recently, the problem of decay of coherent and squeezed packets confined initially 
in a deep potential well of a finite depth was considered by several authors 0, 0, The 
contribution of dissipation was studied, e.g., in |^, ||, |^, |^, ^. However, in all those papers 
the infiuence of the wave packet shape on the transition or escape rates was analyzed in the 
framework of numerical calculations only. 

The aim of the present article is to present simple analytical formulas for the decay rate in 
the case when a wave packet is localized initially near to the bottom of a deep potential well 

V{x) = -muj^ (x^ - Sx") , u>2. (1) 



The special cases of this potential were considered in 0, ^ ||, ||, (z/ = 3) and |TT| (z/ = 4). 



More precisely, we assume that the potential is close to the harmonic one for relatively small 
values of |a;|, while it goes to +oo when x — oo, so that we have a single barrier at a; > 0. 
Besides, it is implied that the potential is given by Eq. provided x"^ — 5x^ > 0, whereas at 
large values of x (to the right of the barrier) it tends to some constant value. 

If the initial state were described by means of the diagonal density matrix p = J2 Pn|^)(^|; 
then the total decay rate 7 would be a sum of partial rates 7^ taken with proper weights 



10, 12, la 



7 = IIPn7n- (2) 

n 

For pure superpositions of many wave functions with different energies En the situation, in 
principle, may be more complicated due to the possibility of quantum interference effects. 
However, under certain conditions Eq. (0) can be applied to the pure states, as well. Suppose 
that we have a single high (unidimensional) barrier from the right, so that the motion can be 
considered as free at x > L. Let us designate the probability of discovering the particle in the 
well and under the barrier as wl = S^^o dx . Then an immediate consequence of the 

quantum continuity equation is the relation wl = —j{L), where j = {iTi/2m) {ipVip* — ip*V'ip) 
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is the usual current density. If the wave function of the packet has the form 

^ = ^c„^„e-^^"*/^ (3) 

n 

(it holds for decaying states, as well, provided that t <^ 7"^), then 



3 ^ ^ I Cji I in ~l~ ^ ^ jmn 



i{E„-E^)t/h 



Since the second sum consists of a large number of rapidly oscillating terms with different 
phases and frequencies, it turns practically into zero for \En — Em\t/h ^ 1, and we get Eq. (Q) 
with pn = |c„p, jn = 7n, and 7 = j, the overbar meaning the averaging over fast oscillations. 
This result holds under the condition h'-f <^ \En — E^l, where suffices n,m correspond to all 
the coefficients c„ that yield significant contributions to the expansion (^). In the special 
case of potential (|1|) with u = 3, Eq. (0) was actually derived in in the framework of the 
quasiclassical method proposed in In the present paper we pursue two main goals: i) to 
find an analytical expression for the partial decay rate 7„ in the potential (|l]), ii) to calculate 
the sum (j^) explicitly for different physically interesting initial wave packets. 

2 Partial decay rates from a parabolic well 



To calculate the partial tunneling rates we use the standard quasiclassical formula P, p!3| 
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7„ = ^exp|--y2^^ ' ^JV{x) - EndxY (4) 



It is convenient to rewrite it in the form 



where Vq is the barrier height: 

2 /2\2/(''"^) 



zu \v J 

Function Fy{t) is given by the integral 

F,{t) = dx[x^ -x" -tf\ (6) 

Ja 
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where a <h are positive solutions to the equation 



(7) 



It is known that for z/ = 3 and v = A the integral (^) can be expressed in terms of the 
complete elliptic integrals of the first and second kind [l^ 

dx 



K(z) 



1 /I — z'^x'^ 



dx. 



'0 ^(1 -x2) (1 -2%2)' '''' V 

At z/ = 4 the roots of Eq. (|^ can be found explicitly, and formula (3.155.9) of yields 



F.{t) = 1 (i + ey"'' (i + e) e (yr^) - ( yr^ 



-3/2 



1 - 



At = 3 we get 



15 



-5/4 



2 1-f -fe^ E K/i_e 1+f K L/l-e 



(9) 



where = (^(j _ — c), and c is the negative root of Eq. (0). Another expression for 7^^) 
was given in 

Assuming that coefficient 6 in the potential energy (|^) is sufficiently small, we have Q = 
VQ/{huj) ^ 1. Then muxl/h = 2vQ/{v — 2) ^ 1, where is the position of the maximum of 
the potential energy. This means that the energies of the low levels practically coincide with 
the harmonic oscillator energy En = huj{n + For these levels En ^ Vq, so we need the 
expansions of the exact expressions and (|^) at t ^ 1. For z/ = 4, the known asymptotics of 



the complete elliptic integrals |14| results in the formula 



(10) 
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For u = 3 the roots of the cubic equation (0) read (to within an accuracy of the order of t^^"^) 



-v^(i-^), b = i-t, e = 2V~t{i-vt). 



However, the expansion of F^it) does not contain odd powers of y^: 

15 4 V64ey 



Both the expressions, ( pUj ) and (|TT]), contain the same term (t/4) ln(t/e). This coincidence is 
not accidental: the leading term in the expansion of the integral (Q) at t — equals (t/4) ln(t/e) 
for any u > 2. Indeed, the contribution to this integral of the domain near the left turning 
point can be represented as 

Fieft = I - t dx^ , 

where A is some finite number. Thus -F/e/t = (V^) sinh(2z) — tz/2 + ■ ■ -, where A = ^/t cosh(2;). 
In the limit of A^/t ^ 1 we get 



left 



A^I2 + (t/4)(lnt - 1) - (t/2) ln(2A) + 0{t^lA^) 



In the vicinity of the right turning point we put x = 1 — e and write 



right 



Jt/(u-2) ^ 3 \1 -2 



As to the integral in the limits from A to B, it has an obvious power expansion with respect 
to t. Consequently, the following expansion holds: 



Mt) = + (t/4) ln(t/e) - fi'k + O {t' \nt) 



where coefficients fj)'^ and f^^^ depend on the concrete value of the exponent v. Then the 
partial tunneling rate reads 



7i^^ = f exp(-^/r + (n + l/2) 

ZTT A, 



Inf '-^ ]+AfW 



Now let us notice that the factor [n + l/2)[ln(r2 + 1/2) — 1] is the leading term of Stirling's 
asymptotical formula 

ln(n!) = (n + l/2)[ln(n + 1/2) - 1] + (1/2) ln(27r) + 0(l/n), 

which works quite well even at n ~ 1. Consequently, under the restriction 

n2ln(Q)/g<l (13) 

the partial decay rates are given by the Poisson distribution: 

7r=7rff' (14) 



Xv 

Specifically, 
/is = 432, 



fi^ = (2A^)-^exp 4/, 



.(1) 



76 



Xu 

uj\ — exp 
V 27r ^ 



i.{o) 



/i4 



64, 



(3) 

7o 



/216 



TT 



Qexp ( 



(4) 

7o 



/32 



— Qexp — —Q . 



TT 
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Strictly speaking, the right-hand side of Eq. contains some additional preexponential 
factor G {Vo,uj, En). But this factor is a smooth function of energy. This means that G ~ 
Go [1 + Cfi'n/Q)], while the leading exponential term was approximated with an accuracy of the 
order of n^lnQ/Q. Consequently, the influence of the preexponential factor can be neglected 
under the restriction ([T^). Note that our expression for ^y^^ coincides identically with the result 
of p|, where a special attention was paid to the correct calculation of the preexponential term. 



Formula (H) seems to be a universal distribution of the partial decay rates from the low 
energy levels, which holds for any potential of the form V{x) = {muj'^ /2)[x'^ — u{x)] with 
\u{x)\ <^ a;^ at a; — > 0, provided that conditions Q ^ 1 and ([T3|) are fulfilled. The concrete 
form of u{x) is responsible for the precise value of the coefficient fi. Since Xu ^ 1? the total 
decay rate turns out very sensitive to the detailes of the distribution pn- 



3 Decay of a slightly deformed ground state 

The simplest example of the initial wave packet corresponds to a coherent state, i.e., an eigen- 
state of the operator a = {mux + ip) / \/2mujTi. Then 

p„(«) = ^e-H, |a|2 = n, (15) 
n! 

and Eqs. (0), ([T^) and (|15]) result in the formula (we drop the suffix p) 

Icoh = 7oexp Jo (2|a|^) = 7oexp (-n) Jq (2^/%^^) , (16) 

Io{z) being the modified Bessel function. To evaluate the domain of its validity, we notice that 
the maximal contribution to the sum (^ is given by the terms with n^ax ~ v'x^, and the 
width of the effective distribution 7nPn is obviously less than ^JUmax- Thus the requirement 
(|13D results in the inequality |ap = n <^ l/lnQ, so that the term exp {—n) can be omitted. 
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However, the condition n <^ 1/ InQ does not exclude the possibihty of i/xn ^ 1. In this case 
we have 

-1/2 



7coh = 7o (4vrv^j ' exp [2^Jxn) > 7o- (17) 
The most general pure Gaussian state (which is called frequently in the current literature 



as a "squeezed state": see [0 and references therein) can be considered as an eigenstate of a 
linear combination of the operators a and a^: 



h\i5uv) = i5\i5uv) , b = ua + va\ - = 1 (18) 

(for simplicity, we confine ourselves to the case of linear uniform transformations). The corre- 



sponding level population distribution reads [|T^, 

1 



pn = I {n\(3uv) 



\u\n\ 



2u 



exp 



\I3\' + Re[(3 



u 



/9 



uv , 



(19) 



where Hn{x) is the Hermite polynomial. To calculate the sum (|^) we need a formula for 
J2 g"'Hn{x)Hn{x*)/{n\)'^. It can be easily found, if one takes the known generating function of 
the Hermite polynomials 

—Hn{x) = exp (2xz - z^] , 

n=0 ^• 

multiplies both sides by the complex conjugated functions, puts z = ^exp(i(p), and integrates 
the product over ip. The result is 



51 r~n2 exp [4\x\y^ cos{!f + '^) -2gcos{2if)] dip, 



X = \x\e 



In this way we get the expression 

7o 



Isq 



27rm JO 



p2n 




d 


/ exp 


2 




Jo 




u 



COs(2(y9) 



dtp, 



if) = arj 



uv 



(20) 



(we neglect the contribution of the terms proportional to in the argument of the exponential 
function, since it is very small under the conditions x ^ 1 ^ind <^ 1, which ensure the validity 



of Eq. (20)). At i; = we arrive again at Eq. (p!6D . 

For a slightly squeezed state with x|t>| <^ |/3|v^ (this inequality implies \v\ <^ 1, so that \u\ 



can be replaced by unity) the integral in (^) can be calculated with the aid of the steepest 
descent method, provided that ^ 1- The integrand assumes its maximal value at 

(p = —ip. Thus we get 

75g(/3, v) = -fcohiP) exp [-xl^^l cos(2^/')] . 



Consequently, the decay rate of a squeezed packet may be both greater and less than the 
decay rate of the coherent packet with the same mean energy, depending on the value of the 
phase difference ■?/'. (This qualitative result was obtained in |15|, 0] for z/ = 3, although the 
quantitative estimations were not quite correct, since the importance of logarithmic terms in 
the expansion of F-^{t) was underestimated.) To elucidate the situation, we take into account 
the formulas for the average number of quanta and its variance in the squeezed state |T^ 



n = \v\ 



(1 + 



[n] 



2n 



2\u\ 



2Re P' — 

u 



They can be simplified significantly if|f|<^l,|M|~l: 

|/3|^[1 -2|t;|cos(2V^)], 



Finally, we get 



Isq = 70 



-1/4 



exp 



2n~\P\\ 



S = {an- n) /n, 



(21) 



where S is the known Mandel's parameter characterizing the type of photon statistics. In this 
case the super-Poissonian statistics enhances the tunneling rate, while the sub-Poissonian one 
suppresses it. A different dependence of the tunneling rate on the degree of squeezing was found 
in Ref. [Q. But its authors performed the numerical calculations in a quite different domain of 
parameters: Q = 10 and \(3\ > 1, where our approach cannot be applied. 

In the case of a squeezed vacuum {(3 = 0) the right-hand side of Eq. ( POD coincides with the 
known integral representation of the modified Bessel function. Then 



'yv = 7o/o(xl^^l) = loHxVn), 



(22) 



provided that |fp = n <^ {QlnQ)^^ (for this reason we put \u\ = 1). If xl^l ^ 1 a-^d 
l/^l I'^lv^' then one can use again the steepest descent method. Now we have two extremal 
points: ip = tt/2 and (/? = 37r/2, so 



7i;(/?) = 7o (SttxI^^I) exp(x|t;|) cosh(2|/3|^sin^). 
In this case we have 7^(/3) > 7,;, although Mandel's parameter 



(23) 



S 



|t;|2-2|/3|>|cos(2V^)] + 



may be both positive and negative, in spite of the requirement \(3\ <^ I'^^ly^- 

The level populations in a Gaussian mixed state with zero mean values of the quadratures 
are expressed in terms of the Legendre polynomials |]18|, |T9|, pO|] : 

n/2 / 

Pn = 2{Ad + 2T + 1] 



-1/2 



Ad + l- 2T 



Ad+1 + 2T 



P. 



Ad - 1 



{Ad + ly - AT^ 



1/2 



(24) 



Parameters d and T are related to the "degree of mixing" of the quantum state and the mean 
quantum number: 



d-' 



T = l + 2n, T> 2Vd. 



(25) 



In this case the sum in the right-hand side of Eq. (^ is reduced to the known generating function 
of the Legendre polynomials (see, e.g., Eq. 10.10(40) from I^Tf). Since we are restricted with 
the inequality n <^ 1, it is convenient to introduce a small parameter e < n <^ 1 according to 
the relations Tr (p^) ^1 — 2e,d^l/A + e. Then we get a simple formula 



Jgaussix, n, e) = 7o exp(x£:)/o [xVn - ej . (26) 

At e = it coincides with (p2D . In the thermal state we have e = n, pn = n"'/ (1 + ^)"^^, and 
Itherm = 7oexp(xn). The last expression holds provided that n <^ (QlnQ)^^/'^ . Note that this 
restriction does not forbid the inequality ^ 1- With the same value of n, at ^ 1 the 
tunneling rate from the squeezed (pure) vacuum state turns out much greater than that from 
the thermal one. However, the thermal state decays faster than the coherent one under the 
same conditions. These examples show that the decay rates are very sensitive to the details of 
the energy distribution in the wave packet, so it is difficult to find a general law. 

An example of a Gaussian packet with nonzero means of the quadratures is the mixture of 
the coherent and thermal states |TB|, |22[ , when the "shifted Planck distribution function" is 
expressed in terms of the Laguerre polynomials: 



n 



Pn 



th 



(1 + ntu)' 



n+l 



exp 







Ln ( 


1 + nth_ 



nth{l + nth] 



(27) 



In this case sum (|]) can be calculated exactly with the aid of Eq. 10.12(18) from |2T[]. For 
nth 1 and |ap -C 1 the total decay rate equals the product of the coherent and thermal 
decay rates: 

lshift{W\,nth) = 7oexp (xnth) h (2|tt|A/x) = Ithermlcoh- (28) 
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An example of a nonGaussian wave packet is the even coherent state introduced in ||23[| , 
\a;+) = {2 [1 + exp {-2\a\^)] } (| a)+ \ -a)) , 
with the quantum distribution function 

I I An 

P2n — (^2?7,)! cosh |aP ' /^2n+i — u, 'i — |«| tann|tt| . 

In this case the total decay rate is proportional to the sum of the usual and the modified Bessel 
functions of the argument 2|Q;|y^, provided that lap <^ (/XylnQ)"^. Then n = \a\^, and 
cr„ = 2n at lal <^ 1. Therefore 



Jo (2ni/^v^) + Jo {2n'/'Vx^ 



7+ = 



We see that the even coherent state is less stable with respect to tunneling than the Glauber 
coherent state with the same value of n -C 1. For all distributions, 7 ~ 70 if 1- 



4 Decay of slightly deformed excited states 

An odd coherent state 

I a; -) = {2 [1 - exp (-2|ap)] }"^^^ (| a)- \ -a)) , 

|^|4n+2 

^^"'^= (2n + l)!sinh|ar = °' ^^"^ = I""' ^^^^ 

is an example of the deformed first excited oscillator state at |a| <^ 1. Its decay rate equals 



7- 



7i 



[Jo i2\a\^) - Jo i2\a\^)] , 71 = 7oX- 



2x|a|2 

The limitations on |a| are the same as above, but now = 3 (n — 1). 



Another example is the odd squeezed state [^: 
\ z; —) = exp 



a' 



1), 



(,) _(2n+l)! 

P2n+1 — 



(nl 



2n 



p2 = 0. 



The expression for the total decay rate is similar to formula (p2D for the squeezed vacuum state: 

72 = 7i^o(xk|), 
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provided that \z\'^ = (n - l)/3 < {QXnQY^. 

It is not difficult to perform the calculations also for two families of deformed m-quantum 
states. The ffist one corresponds to the photon-added coherent states (PACS) ||24|| : 



a,m) = [a}j | a), 



IU|2{"-"*) 



n\\a 



m\[{n — m)!]2 ' 



n > m 



(we assume that \a\ <^ 1). Then the total decay rate is given by the formula similar to (jig) , 



Tpacs = 7m/o (2|a I v^x) , |ap<l/lnQ, 

but with another meaning of parameter |a|, since now n — m = (m + l)|a;p = cr„. 

The second family consists of the displaced number states |25, 26, 27, 28, 2£] | m,a) = 
exp (^aa^ — a*a^ \ m), whose quantum distribution function is expressed in terms of the asso- 
ciated Laguerre polynomials: 



pr" = — r 

ml 



In this case, using the identity 10.12(19) from |21[ 



X 



[i + zT 



n=0 



and applying the same approach that led to Eq. (^OD, one can obtain the formula 



n=0 



X] 



27r Jo 



2tt 



dip exp {—2xy/y cosip) {1 + y + 2y/y cosip)"^ . 



Consequently, the total decay rate can be expressed as some combination of the modified Bessel 
functions of 2|a|y^ with different integer indices. However, in the most interesting case, when 
2|tt|i/x ^ li the steepest descent method leads to a simple formula 



7m(47r|a|yx) exp(2|a|yx), 



sp 



\a\ 



n — m <^ 1/ InQ, an = {2n + l)\a\ 



5 Conclusion 

Two new results seem to be the most important. Firstly, we have found the universal Poisson 
distribution of the partial decay rates from the energy eigenstates in the parabolic potential 
well for a wide class of potential barriers. Secondly, we have demonstrated that the tunneling 
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decay rates are very sensitive to the shape of the wave packet. In particular, if one has initially 
not an exact energy eigenstate, but a combination (pure or mixed) of the states with different 
energies, then the decay rates may be quite different, even when the average energy, coordinate 
and momentum variances, etc., are almost the same. This fact may be important for the 
analysis of various phenomena related to the tunnel effect, when the initial state is not known 
absolutely exactly. 
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